Self-energy values for P states in hydrogen and low-Z hydrogenlike ions 
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We describe a nonperturbative (in Za) numerical evaluation of the one-photon electron self energy 
for 3Pi/2, 3P3/2, 4Pi/2 and 4P 3 / 2 states in hydrogenlike atomic systems with charge numbers Z = 1 
to 5. The numerical results are found to be in agreement with known terms in the expansion of the 
self energy in powers of Za and lead to improved theoretical predictions for the self-energy shift of 
these states. 
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In this brief report, we consider the one-loop self- 
energy shift for excited P states in hydrogenlike systems. 
High-accuracy numerical calculations of this effect are 
notoriously difficult especially in the region of low nuclear 
charge numbers Z, because the renormalization entails a 
loss of up to nine decimal figures in the numerical calcu- 
lation (the bound-electron self-energy is a residual effect 
which corresponds to the difference of the divergent self- 
energy of the bound electron, minus the corresponding, 
also divergent, shift for a free electron). 

Furthermore, the calculations are needed for a 
self-consistent determination of fundamental constants, 
which relies on the experimental/theoretical analysis of 
a number of hydrogenic transitions. 

One may obtain rather accurate values for the so-called 
self-energy remainder functions by an interpolation Q] of 
results for high and low nuclear charge numbers, which 
rely on two different approaches: (i) direct numerical cal- 
culations at high nuclear charge numbers (see e.g. no 
Za-expansion), and (ii) analytic calculations at low nu- 
clear charge numbers, employing the Za-expansion (see 
e.g. 0). In this case, the values obtained for low nuclear 
charge depend on the interpolation method used, as well 
as (of course) on the reliability of both the numerical cal- 
culations at high Z and the analytic calculations for low 
nuclear charge number. 

Here, we follow a third approach and calculate the 
self-energy without Za-expansion, at low Z, using a 
method described previously in jj. 0- El Essentially, 
this method relies on an adequate formulation of the 
physical problem, by which divergent terms are suitably 
identified and calculated separately using semi-analytic 
approaches, and on the use of efficient numerical meth- 
ods for the high-accuracy calculation of the Green func- 
tion of the relativistic electron and for the evaluation of 
slowly convergent sums of intermediate angular momenta 
(it might be useful to mention the keyword "convergence 
acceleration" in that context 0). On modern computer 
processors, it is not even necessary to parallelize the cal- 
culation; that latter approach had previously been em- 
ployed in 0, • 

All calculations reported here are carried out in the 
no-recoil limit, which corresponds to an infinitely heavy 
nucleus. It might be instructive to recall the following 
subtlety, which is well known, and to include a slight 



detour: At the current level of accuracy and especially at 
low Z ', the reduced-mass dependence of the self energy 
should be included for a comparison of the effect with 
experiment. This can be done a posteriori by considering 
the semi-analytic expansion in Eq. J2J) below, using the 
formulas given in Eq. (2.5b) of Ref. @], which indicate the 
reduced-mass dependence of the coefficients, and then the 
self-energy remainder values given in Tables^] — HVI below 
in this brief report. 

Returning to the discussion of our calculation, we write 
the (real part of the) energy shift Ai?sEdue to the elec- 
tron self-energy radiative correction as [g] 

. _ a (Za) 4 m e c 2 . * 

AE SE - - ^ o F(nLj , Za) , (1) 
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where F is a dimensionless quantity (m c is the elec- 
tron mass, a is the fine-structure constant, and c is the 
speed of light in vacuum). In writing the expression 
F(nLj, Za), we follow the usual spectroscopic notation 
for the quantum numbers of the single electron in a hy- 
drogenlike ion; namely, we denote the principal quantum 
number by n, the orbital angular momentum by L and 
the total electron angular momentum by j. 

TABLE I: Numerical results for the scaled self-energy 
function F (SP1/2, Za) and the self-energy remainder 
function Gse, as defined in Eq. |T). in the regime of 
low nuclear charge numbers Z. See Fig. Q 



z 


F(3P 1/2 ,Za) 


G SE (3P 1/2 ,Za) 


1 


-0.115 459 16(5) 


-1.118 5(9) 


2 


-0.114 787 32(5) 


-1.089 7(2) 


3 


-0.113 831 60(5) 


-1.062 5(1) 


4 


-0.112 644 39(5) 


-1.036 30(6) 


5 


-0.111 258 78(5) 


-1.010 84(4) 



The leading terms in the semi-analytic expansion of 
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TABLE II: Same as Table □ for the 3P 3/2 state. See 
also Fig. Q 





r (3*3/2, A<*i 




1 


0.134 414 38(5) 


—0.577 5(9) 


2 


0.134 792 58(5) 


-0.559 3(2) 


3 


0.135 332 45(5) 


-0.542 1(1) 


4 


0.136 006 06(5) 


-0.526 06(6) 


5 


0.136 795 05(5) 


-0.510 76(4) 


TABLE III: Same as Tables Q] and [H] for the 
state. See also Fig.|3 


Z 


F{AP 1/2 ,Za) 


G S e(4P 1/2 ,Zq) 


1 


-0.110 425 6(1) 


-1.164(2) 


2 


-0.109 720 3(1) 


-1.134 1(5) 


3 


-0.108 717 9(1) 


-1.105 5(2) 


4 


-0.107 471 8(1) 


-1.078 1(1) 


5 


-0.106 016 8(1) 


-1.051 20(8) 


F(nPj, Za) about Za — read 




F(nP , 


Za) = A i0 {nPj) + 


(Za) 2 



1/2 



G SE (nP 17 Za)] . (2) 



The A coefficients have two indices, the first of which de- 
notes the power of Za [including those powers explicitly 
shown in Eq. (J)], while the second index denotes the 
power of the logarithm ln(Za)~ 2 . The evaluation of the 



TABLE IV: Numerical results for the scaled self- 
energy function F and the self-energy remainder func- 
tion Gse for the 4P 3 / 2 state. See also Fig. 



z 


F(4P 3/2 ,Za) 


G SE (4P 3/2 ,Za) 


1 


0.139 440 2(1) 


-0.609(2) 


2 


0.139 832 7(1) 


-0.590 4(5) 


3 


0.140 392 7(1) 


-0.572 9(2) 


4 


0.141 091 6(1) 


-0.555 5(1) 


5 


0.141 909 8(1) 


-0.539 41(8) 



coefficient 
has been historically problematic. 



lim Gse(iiS 1 /2, Za) 

Za— >0 



(3) 



c 

N 



(3 



-0.95 I 

-1.0 — 

-1.05 — 
-1.1 — 
1.15 — 
-1.2 — 



2 3 

z 



N 



Q. 

CO 



<3 



-0.6 



2 3 

z 



FIG. 1: Comparison of the self-energy remainder for 3Pi/ 2 
and 3P 3 / 2 states, as listed in Tables HI and ITT1 to their low-Z 
limit which is the A@q coefficient [see Eqs. and ©]. Here, 
Z is the nuclear charge, and the self-energy remainder Gse 
is a dimensionless quantity. The nuclear charge number is 
denoted by Z. 



We now list the analytic coefficients and the Bcthe log- 
arithms relevant to the atomic states under investigation, 
referring the reader to [!j for a more detailed discussion 
and further references, 

1 4 

Ai (nP 1/2 ) = -- - - \nk (nP) , (4a) 

A i0 (nP 3/2 ) = -L _ ! lnfc (nP) . (4b) 

Numerical values for the Bethe logarithms lnfco(nP) are 
well known [l(3,[ll|. 

The A 6 i-coemcients for the states under investigation 
read 



268 

A 6 i(3P 1/2 ) = — , A 61 (3P 3/2 ) 
499 

A 6 i(4P 1/2 ) = — , A 61 (4P 3/2 ) 



148 
405 
137 
360 



(5a) 
(5b) 



The Arc coefficients have been evaluated in 9], and more 



3 



recently in Q to an increased accuracy, 

^6o(3Pi/ 2 ) = -1.148189 956(1), (6a) 

^6o(3-P 3/2 ) = -0.597 569 388(1), (6b) 

A 6 o(4Pi/ 2 ) = -1.195 688142(1), (6c) 

A 60 (4P 3/ 2) = -0.630 945 796(1). (6d) 

Note that the result for 3Pi/ 2 had been given inaccurately 
as -1.14768(1) in 0. 
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FIG. 2: The analog of Fig.dfor 4P 1/2 and 4P 3/2 states. 



Our calculation of the nonperturbative (in Za) elec- 



tron self-energy for the 3Pj state (see Tables[l]aiid|nJ has 
a numerical uncertainty of 2 Hz in atomic hydrogen. For 
the APj states, the numerical uncertainty is l.3xZ 4 Hz 
(see Table ITTT1 and II V|) . In the non-recoil limit, our result 
for 4:P 1 / 2 , Z = 1, corresponds to a self-energy shift of 

A£ SE (4P 1/2 , Z = 1) = -1 404.239(1) kHz , (7) 

which is in agreement with the result — 1 404.240(2)kHz 
obtained in |l2( via an interpolation of the \ow-Z ana- 
lytic results and high-Z numerical data, confirming (in 
this particular case) the validity of the interpolation 
procedure used for various excited hydrogenic states in 
the latest adjustment of the fundamental physical con- 
stants Indeed, all entries for the self-energy re- 
mainder function GW hi Tables III — IIVI are in agreement 
with those used in (lii . IT^j for the determination of the 
fundamental constants, and for the precise calculation 
of hydrogenic energy levels using the method of least 
squares. Our all-order evaluation eliminates any uncer- 
tainty due to the unknown higher-order analytic terms 
that contribute to the bound electron self-energy of 3P 
and 4P states [see Eq. (j2J] and improves our knowledge 
of the spectrum of hydrogenlike atoms (e.g. atomic hy- 
drogen, He + ). Furthermore, the numerical data for the 
self-energy remainders check the validity of the highly in- 
volved analytic approach that led to the evaluation of the 
Aeo-coefncients as listed in Eq. (JHJ). 
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